Comparison of Source Images for protons, 7r 's and A's in 6 AGeV Au+Au collisions 
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Source images are extracted from two-particle correlations constructed from strange and non- 
strange hadrons produced in 6 AGeV Au + Au collisions. Very different source images result from 
pp vs pA vs n~ir~ correlations. These observations suggest important differences in the space-time 
emission histories for protons, pions and neutral strange baryons produced in the same events. 



Rclativistic heavy ion collisions of 1-10 AGeV pro- 
duce a fireball of nuclear matter with extremely high 
baryon and energy density 0. The dynamical evolution 
of this fireball is driven by such fundamental properties 
as the nuclear Equation of State (EOS) and possibly 
by a phase transition, e.g., to a Quark Gluon Plasma 
(QGP) @, @, H H @ Two-particle correlation studies, 
for various particle species, provide an important probe 
of the space-time extent of this fireball 0, H, 0, E3- ^- e " 
cent model calculations suggest that the time scale for 
freeze-out of strange and multi-strange particles ma y b e 
much shorter than that for non-strange particles |lllll2l |. 
implying a much smaller space-time emission zone for 
strange particles. If this is indeed the case, then cor- 
relation studies involving strange particles may serve as 
important "signposts" for dynamical back-tracking into 
the very early stage of the collision where large energy 
densities are achieved . 



In this letter we compare the source properties for pro- 
tons, 7r _ 's and A hyperons extracted from pp, ir~ir~ and 
pA correlation functions, as produced in 6 AGeV Au+Au 
collisions. These data are unique in that they constitute 
the first measurement of pA correlations. If A hyperons 
are in fact emitted from a source with a smaller space- 
time extent, then between this and ir~ source broadening 
from resonance decays, one might naively expect an or- 
dering of two-particle source sizes: R p \ < R PP < Rttw 
On the other hand, at these energies the 3-dimensional 
7r~ radii exhibit mj- scaling |?1 and this should mani- 
fest itself in the angle-averaged n~ sources. Furthermore, 
since hit scaling can be ascribed to position-momentum 
correlations in the particle emission function |l5j . one 
might expect similar effects in the pp and pA sources. As 
we will show, neither an interpretation based solely on 
position-momentum correlations nor on naive geometri- 
cal arguments can fully account for the relative size of 
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the three sources. 

Traditional correlation analyses rely on the weakness 
of Final-State Interactions (FSI). With this, one may 
correct for FSI leaving a correlation that is the Fourier 
Transform of the two-particle source function. This ap- 
proximation is mostly valid for pions; for massive and/or 
strongly interacting particles, such as protons or A hyper- 
ons, this approximation breaks down. Recently, Brown 
and Danielewicz have presented an i mag ing technique for 
analyzing two-particle correlations |16| . The technique 
actually uses the FSI, encoded in the form of the parti- 
cles' final state wavefunction, to extract the two-particle 
source function directly The imaging technique has 
been used to address only a few data sets; Panitkin et 
al.\v^ have shown that this approach gives source radii 
similar to the conventional HBT approach (under the as- 
sumption of a Gaussian source) for tt~it~ pairs emitted 
in central collisions for 2-8 AGeV Au+Au. By contrast, 
Verde et al. 0] have shown very different results for pp 
pairs from 75 AMeV 14 N+ 197 Au. For our purpose, the 
imaging technique's key feature is that one can easily 
compare source functions across different species and as- 
sess the different space-time scales relevant for each par- 
ticle pair - a feature that has been mentioned 0], but 
never seriously utilized. 

We use the imaging technique on pp, tt~it~, and 
pA pairs to extract and compare the source proper- 
ties for protons, 7r - 's and A hyperons. The measure- 
ments have been performed with the E895 detector at the 
Brookhaven Alternating Gradient Synchrotron. Here, we 
concentrate on the construction of the pA and pp correla- 
tions and the results of the imaging analysis. Details on 
the detector and its setup have been reported elsewhere 

We reconstructed the A's from the daughters of their 
charged particle decay, A — ► p + ir~ (Branchin g Ra tio ~ 
64 %), following the procedure outlined in Refs.lljj. Fig- 
urenshows the invariant mass spectrum for A's obtained 
in semi-central (upper 23% of total inelastic cross-section, 
b <7 fm) 6 AGeV Au + Au collisions. For the pA corre- 
lation analysis, an enriched sample (~ 80%) of A's with 
1.11 < M A < 1.122 GcV was used. 

Figures I3a),(b),(c) show the correlation functions, 
C(q), obtained by taking the ratio of foreground to back- 
ground distributions in relative-momentum for Ap, pp 
and n~ir~ pairs respectively. Here, q = \\/—{pi — P2) 2 
is half of the relative momentum between the two parti- 
cles in the pair cm. frame. We applied no explicit gates 
on transverse momentum and rapidity except those im- 
plicit in the tracking acceptance. The mean transverse 
momenta (Pt) for particles from pairs with q < 50McV 
were 0.3GeV and 0.12GeV for pp and ir~ir~ pairs and 
the mean rapidities (y) were -0.33 and 0.1 . Ap pairs 
with q < 50MeV have (P t ) = 0.46GeV and (y)=-0.18 for 
A's and (P t ) = 0.38GeV and (j/)=-0.18 for protons. 

We constructed the numerator (or foreground) distri- 
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FIG. 1: A invariant mass spectrum from semi-central 6AGeV 
Au+Au collisions. The hatched area depicts the A mass gate 
used to select the A's for Ap correlation analysis. 
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FIG. 2: Raw smeared (filled symbols) and corrected un- 
smeared (open) correlation functions, C(q), for (a) Ap, (b) 
pp and (c) n~ir~ from 6AGeV Au+Au collisions (b<7fm). 



bution from pairs of particles from the same event, and 
obtained the denominator (or background) distribution 
by pairing particles from different events. We used a 
track- merging filter similar to that outlined in Ref.Q 
to eliminate possible distortions resulting from track- 
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merging effects in the TPC. For each correlation function, 
we used the accepted range of particle multiplicities to 
specify impact parameters b < 7 fm. This range was cho- 
sen to optimize the statistical significance for pA pairs. 

We utilized approximately 100,000 A's (80% purity) 
to yield about 31,000 Ap pairs in the foreground distri- 
bution, with q < 100 MeV/c, after applying the track- 
merging cut for the 6 AGeV data. We have corrected the 
Ap correlation function for: (a) the combinatoric back- 
ground (~ 20%) included in the A sample, (b) feed-down 
due to the electromagnetic decay of the E° (estimated 
to be ~ 25% from RQMD calculations) and (c) smearing 
due to momentum resolution of the detector. 

We corrected the sources for momentum resolution us- 
ing two independent methods, both leading to consis- 
tent results. In the first method, we left the data un- 
corrected and modified the kernel used in the imaging 
analysis to include the smearing effect, assuming an av- 
erage 8p for each (vt) and (y). This technique will be 
detailed elsewhere |2 ij . In the second method, we first 
corrected the correlation functions via an iterative proce- 
dure and imaged with an unsmeared kernel. We explain 
this procedure in the next paragraph. We determined 
the momentum resolution for ir~ and p from GEANT 
simulations giving an average value of 2.0% and 3.1% for 
each component of the tt~ and proton momentum reso- 
lution, respectively. For the A momentum resolution, we 
extracted an average value of 4.0% for each component 
of the resolution using the width of the S~ mass peak 
(S _ decays to Att~) and the ir~ momentum resolution. 

Our iterative procedure starts with model calculations 
(Gaussian sources) for the particle momenta and their 
(unsmeared) correlations, C„(q). We then use a Monte 
Carlo method to smear their momenta and use these 
to calculate the smeared correlation C s (q). The ratio 
C s (q)/C u (q) is used to make a first-try correction to the 
raw observed data correlation C(q). This gives a first it- 
eration unsmeared correlation C' u (q). The latter is then 
smeared in the second iteration (using pairs of particle 
from mixed data events) to give a second smeared corre- 
lation, Cg(q). Typically, the comparison between the raw 
observed C(q) in the data and the second smeared corre- 
lation function led to a reduced \ 2 — !• The associated 
function C' u (q) was then taken as the unsmeared corre- 
lation for the data. Figure |5] shows both smeared and 
unsmeared correlation functions. The results presented 
in this paper are from the iterative method. 

pA correlations can lead to residual correlations be- 
tween primary protons and daughter protons from A de- 
cays. Wang has calculated the magnitude of this residual 
effect on the pp correlation j2^| ■ We have determined its 
effect on our pp correlations to be negligible. This is 
due to (1) the intrinsically small residual pp correlation 
(maximum 3%) arising from our observed pA correlation 
and (2) the small fraction of secondary protons resulting 
from A decay (« 6% of the total number of protons) . The 
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FIG. 3: Correlation functions, C(q), for pA, pp and n~n~ 
pairs are shown in panels (a), (b) and (c) respectively. The 
corresponding short-range source functions, S(r), are indi- 
cated in panels (d), (e) and (f). As a consistency test, a 
simulated correlation function (open squares, circles and tri- 
angles) is recalculated from S(r). Hatched bands show the 
zone of one standard deviation. 

same reasoning applies for the expected perturbation of 
pA correlations from AA correlations [23j and of 7r~7r~ 
correlations from 7r~A or ir~ K° correlations [24(. 

Figures. |^a),(b),(c) show correlation functions for Ap, 
pp and 7r~7r~ pairs respectively. We have not corrected 
the pp and tt~t:~ correlation functions for the Coulomb 
interaction since this effect is included in the imaging pro- 
cedure. The two particle correlation and the source func- 
tion are related through the Koonin-Pratt equation [25) : 

C{q) - 1 = 4tt J drr 2 K(q 7 r)S(r). (1) 

This is a linear integral equation that we may invert to 
obtain the source function S(r) using the techniques in 
Ref. 01 . Here, the imaged source function S(r) gives 
the probability of emitting a pair of particles a distance 
r apart in the pair cm. frame. The derived source 
functions are shown in Figs. El (d),(e),(f). As a consis- 
tency check, we have recalculated the correlation func- 
tions from the derived source functions and these are 
shown in Figs. E[a),(b), ( c ) as restored correlations. 

In Eq. Q), the kernel K(q,r) encodes the FSI and is 
given in terms of the final state wavefunction as K(q, r) = 
\ J d(cos(0 qi r))(|$ q (r)| 2 — 1). In this work, we used the 
Coulomb force for the pion and proton pairs. Addition- 
ally, we used the Reid93 nucleon-nucleon force for pro- 
tons [2(J and the phenomenological potential of Bodmer 
and Usmani for Ap pairs (22) ■ Since the Ap potential 
is not very well known, we reanalyzed the pA correla- 
tion using a simplified kernel that depends only on the 
pA effective range and scattering length. We found no 
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significant change in the imaged Ap source. 

The sources in Figs. (d),(e),(f) appear Gaussian, al- 
though we cannot definitively conclude this given the size 
of the error bands on the imaged sources. In principle, 
the source function can be composed of an admixture of 
short and large-range emission sources 0. In practice, 
the shape of the correlation is strongly dominated by 
the short-range source and the less-corrclated pairs from 
any large-range source essentially dilute the strength of 
the observed correlation. Let us then assume a Gaussian 
for the short-distance part of the time integrated emis- 
sion function for particle type i: Di ~ /^cxp (— r 2 /2i? 2 ). 
The fraction of particles emitted from this source is 
< fi < 1. This choice gives us a Gaussian two-particle 
source function for particles i and j: 



^(r)c<A^exp(-r 2 /2(i? l 2 + i?f)) 



(2) 



To ensure that we use the same source parameterization 
for like and unlike pairs, we choose R 2 tj = |(JJf + R 2 ). 
With this choice the emission function radius is exactly 
the two particle source radius if i = j (i.e. Ru = Ri) |28j . 
Here the fraction of pairs \j is related to the fraction of 
particles in each emission function through \j = fifj. In 
all subsequent discussion, we take the R\/2 value (the ra- 
dius at half maximum density) directly from the sources, 
but we convert the source height into an equivalent Gaus- 
sian A. Values for A and R\/2 of the short-range sources 
are indicated in Fig.0 (d),(e),(f). 

Scanning Fig. EId),(e),(f), the pion source is clearly 
the broadest. However, our naive expectation that the 
Ap source would be the smallest source appears wrong. 
In an effort to understand why, we investigate the source 
sizes in more detail. 

The results in Fig. E£c),(f) for 7r~7r~ show that R\/2 
= 8.2±0.6 fm. This radius, corresponding to a centrality 
range of b < 7 fm, is identical to that reported by Pan- 
itkin et al for b < 5 fm|T3,E3- The A value of 0.25±0.03 
indicates that half of the pions arise from a source with 
R\/2 ~ 8.2 fm and the other half from a considerably 
larger source, possibly caused by resonance decays. For 
pA pairs [cf. Figs. EJa),(d)] the correlations are domi- 
nated by a source of smaller size (R1/2 = 4.8±0.8 fm) 
comprising 60% ± 20% of the pA pairs. 

The contrast in size with the pp source is pronounced 
(cf. Figs. EIb),(c)). The pp correlations are dominated 



by a very compact source Ri 



/2 



3.3±0.2 fm and this 



source comprises only 26% of the pairs or 51% of the to- 
tal protons. Both in 75 AMcV 14 N+ 197 Au Q and 200 
AGeV S+Pb investigators have also found relatively 
small source sizes (R1/2 = 2.9 fm and R1/2 = 3.2 fm, 
respectively) for roughly the same fraction of protons. 
From these observations, we suspect that the there is a 
common feature in nucleus-nucleus collisions that span a 
broad range because (1) strong collective motion focuses 
the source to much smaller radii |30| . and (2) pp corre- 
lations are insensitive to the long distance parts of the 
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FIG. 4: Variation of source half-radius R1/2 (top panels) and 
A parameter (bottom panels) for different phase space regions 
identified by their mean transverse momentum (Pt) and mean 
rapidity (y) values (shown in different symbols) for sources 
from (a) and (b) pp correlations. The open cross in 

figure (a) shows the R1/2 value, calculated from the HBT 
radii from Ref [7]. 



source, as discussed by Wang and Pratt [13|. While we 
cannot separate these two effects in this study (or even 
rule out other more exotic causes) we comment that pi- 
ons are much less focused by collective effects because 
of their much lower mass [3l|. Consequently, we expect 
the pion correlation function to probe a larger portion of 
the emission region and hence be associated with a larger 
source. 

Previous works 0, 0] have indicated that the 3- 
dimensional pion radii exhibit tut scaling attributable to 
collective flow. R1/2 in the pair cm. frame can be com- 
puted from the 3-dimensional source radii i? ut , Rside and 
i?iong in the Longitudinally Comoving Source (LCMS) 
frame via 



R 1/2 « 1.66{/7i?out-Rsidc#i 



ong, 



(3) 



where the Lorentz factor is 7 = mx/m. In the various 
models of m<r scaling of correlation radii, the radii in the 
LCMS frame all have the same form: R = R y/T/mT- 
Here T is the source temperature and the R 0, s for each 
radius depend on the details of the model in question. 
Inserting this dependence into Eq. Q gives 



r 1/2 « BP 1/2 Jrjrt 



11 

T ■ 



(4) 



where 



,cir 



rriTm 2 . This implies that R1/2 for a 



given system scales weakly as m T '' 6 regardless of the 
rriT scaling model as is evident in Fig. 4. 

If one assumes that R®/ 2 and T are the same for both 
7r~ and proton sources then Eq. (^predicts the ratio 
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Rpp/R^-^- = y"m^_/m|* = 0.4 which is consistent 

with the value of the experimental ratio (0.40 ± 0.04). 

In contrast to the pp and tt~tt~ sources, the pA source 
size does not show this scaling. Naively applying 
Eq. Q and keeping a constant temperature and geomet- 
rical source size, one obtains a predicted pA source size 
of 3.1 ± 0.2 fm which is significantly smaller than the ex- 
perimental value. This value cannot be accounted for via 
the proton-A mass difference of ss20%. It could reflect 
differences between the emission time and the extent of 
collective focusing for A's and protons. In fact, the mag- 
nitude of collective flow for A's is known to be smaller 
than that for protons ■ 

As a possibility, let us assume that the underlying pro- 
ton emission function is common to both the pp and 
pA correlations. For this assumption to hold, we must 
assume that any flow induced focusing affects the pp 
and pA sources similarly and that we may concentrate 
on only the short-range proton and A sources. In this 
context, we can use Eq. and the pp source function 
values to extract information about the A emission func- 
tion. We find that / A = \h P /^f\p~p~ = 1.18 ± 0.40 and 

R\ = ^2R 2 Ap - R? pp = 5.93±1.30 fm. Given that f\~l, 
we may argue that all A's are made from this source. 
Such a moderate sized A source could indicate the ge- 
ometrical size of the hot central region of the collision 
zone, where the energy density is high enough for the 
production of the strange quarks that form the observed 
A's. 

In summary, we have measured small- angle correla- 
tions for pA, pp and tt~tt~ pairs produced in 6 AGeV 
Au+Au reactions and have analyzed them by the source 
imaging technique |l6l | . The strong differences in effective 
source radii reflect very different dynamical histories for 
each pair. The pp and tt~~tt~ pairs may reflect a small 
homogeneity length caused by flow focusing the source 
while the A emission functions reflect a spread out par- 
ticipant zone. 
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